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We discuss how to construct models of interacting anyons by generalizing quantum spin 
Hamiltonians to anyonic degrees of freedom. The simplest interactions energetically favor 
pairs of anyons to fuse into the trivial ( "identity" ) channel, similar to the quantum Heisenberg 
model favoring pairs of spins to form spin singlets. We present an introduction to the theory 
of anyons and discuss in detail how basis sets and matrix representations of the interaction 
terms can be obtained, using non-Abelian Fibonacci anyons as example. Besides discussing 
the "golden chain" , a one-dimensional system of anyons with nearest neighbor interactions, 
we also present the derivation of more complicated interaction terms, such as three-anyon 
interactions in the spirit of the Majumdar- Ghosh spin chain, longer range interactions and 
two-leg ladders. We also discuss generalizations to anyons with general non-Abelian SU(2)fc 
statistics. The k — > oo limit of the latter yields ordinary SU(2) spin chains. 

§1. Introduction 

While in classical mechanics the exchange of two identical particles does not 
change the underlying state, quantum mechanics allows for more complex behav- 
ior. In three-dimensional quantum systems the exchange of two identical particles 
may result in a sign-change of the wavefunction which distinguishes fermions from 
bosons. Two-dimensional quantum systems - such as electrons confined between 
layers of semiconductors - can give rise to exotic particle statistics, where the ex- 
change of two identical (quasi) particles can in general be described by either Abelian 
or non-Abelian statistics. In the former, the exchange of two particles gives rise to 
a complex phase e* e , where 9 = 0,7r correspond to the statistics of bosons and 
fermions, respectively, and 9 ^ 0, tt is referred to as the statistics of Abelian anyons. 
The statistics of non-Abelian anyons are described by k x k unitary matrices acting 
on a degenerate ground-state manifold with k > 1. In general, two such unitary 
matrices A, B do not necessarily commute, i.e. AB ^ BA, or in more mathematical 
language, the k x k unitary matrices form a non-Abelian group when k > 1, hence 
the term non-Abelian anyons. 

Anyons appear as emergent quasiparticles in fractional quantum Hall states and 
as excitations in microscopic models of frustrated quantum magnets that harbor 
topological quantum liquids. 1 ) -3 ) While for most quantum Hall states the exchange 
statistics is Abelian, there are quantum Hall states at certain filling fractions, e.g. 
v = | and v = t 2 -, for which non-Abelian quasiparticle statistics have been pro- 
posed, 4 )' 5 ) namely those of so-called Ising anyons 6 ) and Fibonacci anyons 7 ) respec- 
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tively. Non-Abelian anyons have also generated considerable interest in proposals for 
topological quantum computation, 8 ) where braiding of anyons is used to perform the 
unitary transformations of a quantum computation. The simplest anyons with non- 
Abelian braiding statistics that can give rise to universal quantum computational 
are the so-called Fibonacci anyons which we will discuss in detail in this manuscript. 

In the following, we will first give a short introduction to the mathematical 
theory of anyons, and discuss how to (consistently) describe the degenerate manifold 
of a set of (non-interacting) anyons. Having established the basic formalism we will 
then turn to the question of how to model interactions between anyons and explicitly 
construct matrix representations of generalized quantum spin Hamiltonians. We then 
discuss an alternative formulation in terms of non-Abelian SU(2)& anyons. Rounding 
off the manuscript, we shortly review some recent work analyzing the ground-state 
phase diagrams of these Hamiltonians. 



§2. Basic theory 



2.1. Algebraic theory of anyons 

In general terms, we can describe anyons by a mathematical framework called 
tensor category theory. In such a categorical description, anyons are simple objects 
in the corresponding tensor categories, and anyon types are the isomorphism classes 
of anyons. Here we will not delve into this difficult mathematical subject, but focus 
on the theory of Fibonacci anyons, where many simplifications occur. 

2.2. Particle types and fusion rules 

To describe a system of anyons, we list the species of the anyons in the system, 
also called the particle types or topological charges or simply labels (and many other 
names); we also specify the anti-particle type of each particle type. We will list the 
particle types as {xi}™~Q, and use {Xi}™~Q to denote a representative set of anyons, 
where the type of Xj is X{. 

In any anyonic system, we always have a trivial particle type denoted by 1, which 
represents the ground states of the system or the vacuum. In the list of particle types 
above, we assume xq = 1. The trivial particle is its own anti-particle. The anti- 
particle of Xi , denoted as X* , is always of the type of another Xj . If X{ and X* are 
of the same type, we say Xi is self-dual. 

To have a non-trivial anyonic system, we need at least one more particle type 
besides 1. The Fibonacci anyonic system is such an anyonic system with only two 
particle types: the trivial type 1, and the nontrivial type r. Anyons of type r are 
called the Fibonacci anyons. Fibonacci anyons are self-dual: the anti-particle type 
of t is also r. Strictly speaking, we need to distinguish between anyons and their 
types. For Fibonacci anyons, this distinction is unnecessary. Therefore, we will refer 
to r both as an anyon and its type, and no confusion should arise. 

Anyons can be combined in a process called the fusion of anyons, which is similar 



*' Roughly speaking, a universal quantum computer is a general-purpose quantum computer 
which is capable of simulating any program on another quantum computer. 
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to combining two quantum spins to form a new total spin. Repeated fusions of the 
same two anyons do not necessarily result in an anyon of the same type: the resulting 
anyons may be of several different types each with certain probabilities (determined 
by the theory). In this sense we can also think of fusion as a measurement. It follows 
that given two anyons X, Y of type x, y, the particle type of the fusion, denoted as 
X ® Y, is in general not well-defined. 

Given an anyon X, if the fusion of X with any other anyon Y (maybe X itself) 
always produces an anyon of the same type, then X is called an Abelian anyon. If 
neither X nor Y is Abelian, then there will be anyons of more than one type as 
the possible fusion results. When such fusion occurs, we say that the fusion has 
multi-fusion channels. 

Given two anyons X, Y, we formally write the fusion result as X <S> Y = (BiUiXi, 
where Xj are all anyons in a representative set, and n$ are non- negative integers. 
The non- negative integer rtj is called the multiplicity of the occurrence of anyon Xj. 
Multi- fusion channels correspond to Y2i n i > !• Given an anyonic system with anyon 
representative set -pQ}£=o > then we have Xj tg> Xj = ®^~ Q N^-Xk, or equivalently, 
x% (g> Xj = (Bk~ Nj°jXk. The non-negative integers N-*j are called the fusion rules of 
the anyonic system. If Nf j ^ 0, we say the fusion of X{ and Xj to X^ is admissible. 

The trivial particle is Abelian as the fusion of the trivial particle with any other 
particle X does not change the type of X, i.e., 1 ® x = x for any type x. 

For the Fibonacci anyonic system the particle types are denoted as 1 and r, and 
the fusion rules are given by: 

1 T = t 

T (g) 1 = T 
T®T = 10T, 

where the © denotes the two possible fusion channels. 

2.3. Many anyon states and fusion tree basis 

A defining feature of non-Abelian anyons is the existence of multi-fusion chan- 
nels. Suppose we have three r anyons localized in the plane, well-separated, and 
numbered as 1,2,3. We would like to know when all three anyons are brought to- 
gether to fuse, what kinds of anyons will this fusion result in? When anyons 1 and 2 
are combined, we may see 1 or r. If the resulting anyon were 1, then after combin- 
ing with the third r, we would have a r anyon. If the resulting anyon were r, then 
fusion with the third anyon would result in either 1 or r. Hence the fusion result is 
not unique. Moreover, even if we fix the resulting outcome as r, there are still two 
possible fusion paths: the first two r's were fused to 1, then fused with the third r 
to r, or the first two r's were fused to r, then fused with the third r to r. Each such 
fusion path will be recorded by a graphical notation of the fusion tree, see Fig. [TJ 

A fusion path is a labeling of the fusion tree where each edge is labeled by a 
particle type, and the three labels around any trivalent vertex represent a fusion 
admissible by the fusion rules. If not all particles are self-dual, then the edges of 
the fusion tree should be oriented. We always draw anyons to be fused on a straight 
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Fig. 1. Fusion trees of three Fibonacci anyons (top row). 

line, and the fusion tree goes downward. The top edges are labeled by the anyons 
to be fused, and the bottom edge represents the fusion result and is also called the 
total charge of the fused anyons. 

In general, given n r-anyons in the plane localized at certain well separated 
places, we will assume the total charge at the oo boundary is either 1 or r. In 
theory any superposition of 1 or r is possible for the total charge, but it is physically 
reasonable to assume that such superpositions will decohere into a particular anyon 
if left alone. Let us arrange the n anyons on the real axis of the plane, numbered as 
1, 2, • • • , n. When we fuse the anyons 1, 2, • • • , n consecutively, we have a fusion tree 
as below: 



T T T T T T 




The ground-state manifold of a multi-anyon system in the plane even when 
the positions of the anyons are fixed might be degenerate: there is more than one 
ground state (in reality the energy differences between the different ground states 
go to exponentially as the anyon separations go to infinity; we will ignore such 
considerations here, and always assume that anyons are well separated until they 
are brought together for fusion.) Such a degeneracy is in fact necessary for non- 
Abelian statistics to occur. How can we describe a basis for this degenerate ground 
state manifold? 

As we see in the example of three r anyons, there are multi-fusion paths, which 
are represented by labelings of the fusion tree. We claim that these fusion paths 
represent an orthonormal basis of the degenerate ground-state manifold PI 

*' We will not further justify this assertion, but mention that in the conformal field theory 
(CFT) description of fractional quantum Hall liquids the ground states can be described by conformal 
blocks, which form a basis of the modular functor. Conformal blocks are known to be represented 
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The fusion tree basis of a multi-anyon system then leads to a combinatorial 
way to compute the degeneracy: count the number of labelings of the fusion tree or 
equivalently the number of fusion paths. Consider n T-anyons in the plane with total 
charge r, and denote the ground state degeneracy as F n . Simple counting shows that 
Fq = and F\ = 1. Easy induction then gives F n+ \ = F n + F n -%. This is exactly 
the Fibonacci sequence, hence the name of Fibonacci anyons. 

As alluded to above, when two r anyons are fused, 1 and r each occurs with a 
certain probability. This probability is given by the so-called quantum dimension of 
an anyon. Consider the fusion coefficients Nf- of a theory, if we regard the particle 
types Xi as variables and the fusion rules as equations for Xj. Then in a unitary 
theory the solutions di of Xi which are > 1 are the quantum dimensions of the 
anyons of type Xj. di is also the Perron- Frobenius eigenvalue of the matrix Ni whose 

(j, fc)th entry is N^. We also introduce the total quantum order D = yJ^ i df. The 
quantum dimension of the trivial type 1 is always do = 1. In the Fibonacci theory, 
the quantum dimension of r is the golden ratio ip = 1+ 2 V ^ ■ ■ When two r anyons fuse, 
the probability to see 1 is po = A , and the probability to see r is p\ = = ~ . 

2.4. F-matrices and pentagons 

In the discussion of the fusion tree basis above, we fuse the anyons 1, 2, • • • , n 
consecutively one by one from left to right, e.g., n = 3 gives the left fusion tree 
below. We may as well choose any other order for the fusions. For example, in the 
case of three r's with total charge r, we may first fuse the second and third r's, then 
fuse the resulting anyon with the first r. This will lead to the fusion tree on the 
right as shown in Fig. [21 




Fig. 2. (color online) The two fusion trees of three anyons that both result in the same anyon d are 
related by an 11 F -move" . 

Given n anyons with a certain total charge, then each order of the fusions is 
represented by a fusion tree, and the admissible labelings of the respective fusion 
trees each constitute a basis of the multi-anyon system. 

The change from the left fusion tree to the right fusion tree in Fig. [5] is called 
the F-move. Since both fusion tree bases describe the same degenerate ground state 
manifold of 3 anyons with a certain total charge, they should be related by a unitary 
transformation. The associated unitary matrix is called the F-matrix. The F-matrix 
will be denoted as F% bc , where a, b, c are the anyons to be fused, and d is the resulting 
anyon or total charge (Complications from fusion coefficients Nf- > 1 are ignored.) 

For more than three anyons, there will be many more different fusion trees. To 



by labeled fusion trees, which we refer to as fusion paths. 
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have a consistent theory, a priori we need to specify the change of basis matrices 
for any number of anyons in a consistent way: for example as shown in Fig. [3] the 
left-most and right-most fusion trees of four anyons can be related to each other by 
.F-moves in two different sequences of applications of F- moves. 




Fig. 3. (color online) The pentagon relation for the "F-moves". 



Fortunately, a mathematical theorem guarantees that the consistency equations 
for the above fusion trees, called the pentagons, are all the equations that need to be 
satisfied, i.e., all other consistencies are consequences of the pentagons. Note that 
the pentagons are just polynomial equations for the entries of the F-matrices. 

To set up the pentagons, we need to explain the consistency of fusion tree bases 
for any number of anyons. Consider a fusion tree T, and a decomposition of T 
into two sub- fusion trees T\,T% by cutting an edge; the resulting new edge of Ti,T2 
will also be referred to as edge e. The fusion tree basis for T has a corresponding 
decomposition: if Xj's are the particle types of the theory (we assume they are all 
self-dual), for each Xi, we have a fusion tree basis for Ti,T2 with the edge e labeled 
by Xi. Then the fusion tree basis of T is the direct sum over all Xj of the tensor 
product: (the fusion tree basis of Ti) (g) (the fusion tree basis of T2). 

In the pentagons, an F-move is applied to part of the fusion trees in each step. 
The fusion tree decomposes into two pieces: the part where the F-move applies, and 
the remaining part. It follows that the fusion tree basis decomposes as a direct sum 
of two terms: corresponding to 1 and r. 

Given a set of fusion rules Nf- solving the pentagons turns out to be a difficult 
task (even with the help of computers). However, certain normalizations can be 
made to simplify the solutions. If one of the indices of the F-matrix a, b, c is the 
trivial type 1, we may assume F^ ,b ' c = 1. In the Fibonacci theory, we may also 
assume F^ ,b ' c = 1. It follows that the only non-trivial F-matrix is Fr' T,T , which is a 
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2x2 unitary matrix. 

There are many pentagons even for the Fibonacci theory depending on the four 
anyons to be fused and their total charges: a priori 2 5 = 32. It is easy to see that the 
only non-trivial pentagon for F is the one with 5 r's at all outer edges. The pentagon 
is a matrix equation for F extended to a bigger Hilbert space. To write down the 
pentagon, we need to order the fusion tree basis with respect to the decomposition 
above carefully. 

Written explicitly for Fibonacci anyons the pentagon equation reads 

(Fr) d a (fdi = (ffi: (K eT )i (fd: , e-i) 

where the indices a, b, c, d, e label the inner edges of the fusion tree as shown in 
Fig. There are only a few different matrices appearing, of which four are uniquely 
determined by the fusion rules 

prrr = plrr = prlr = prrl = j ^ 

in a basis {1, r} for the labeling on the central edge. The only nontrivial matrix is 
F^ TT . Setting b = c = 1 the pentagon equation simplifies to 

(FD^iFD'AFni, (2-3) 

which combined with the condition that F^ TT is unitary constrains the matrix, up 
to arbitrary phases, to be 

where ip = (y/E + l)/2 is the golden ratio. 

2.5. R-matrix and hexagons 

Given n anyons Y{ in a surface S, well-separated at fixed locations pi, we may 
consider the ground states V(S;pi,Y) of this quantum system. Since an energy 
gap in an anyonic system is always assumed, if two well-separated anyons Yi,Yj 
are exchanged slowly enough, the system will remain in the ground states manifold 
V(S;pi, Yj). If £ V(S;pi, Yi) is the initial ground state, then after the exchange, 
or the braiding of the two anyons Y^, Yj in mathematical parlor, the system will be 
in another ground state |fiq) = Yli^i e i m V{S\Pi->Yi)-> where ej is an othonormal 
basis of the ground states manifold V(S;pi,Y)- When \]Pq) runs over the basis e», 
we obtain a unitary matrix Rij from V(S;Pi,Yi) to itself. In mathematical terms, 
we obtain a representation of the mapping class group of the punctured surface S. 
If S is the disk, the mapping class group is called the braid group. In a nice basis of 
V(S;pi,Y), the braiding matrix Rij becomes diagonal. 

To describe braidings carefully, we introduce some conventions. When we ex- 
change two anyons a, b in the plane, there are two different exchanges which are 
not topologically equivalent: their world lines are given by the following two pic- 
tures, which are called braids mathematically. In our convention time goes upwards. 



8 



S. Trebst, M. Troyer, Z. Wang, and A. Ludwig 



U 

n 



u 
n 



When we exchange two anyons, we will refer to the right process, which is called the 
right-handed braiding. The left process is the inverse, left-handed braiding. 

Now a comment about fusion trees is necessary. In our convention, we draw the 
fusion trees downwards. If we want to interpret a fusion tree as a physical process 
in time, we should also introduce the conjugate operator of the fusion: splitting of 
anyons from one to two. Then as time goes upwards, a fusion tree can be interpreted 
as a splitting of one anyon into many. 

All the braiding matrices can be obtained from the i?-matrices combined with 
.F-matrices. Let Vc ,b be the ground state manifold of two anyons of types a, b with 
total charge c. Let us assume all spaces Vc ,b are one-dimensional, and et' b be its 
fusion tree basis. 

When anyons a and b are braided by R a ,b, the state e"' 6 in Vc' b is changed into 
a state R a f,&c b in V b,a . Since both R a ^et' b and e b ' a are non-zero vectors in a one- 
dimensional Hilbert space V c ' a , they are equal up to a phase, denoted as R b ' a , i.e, 
Ra,b e< c = Rc ,a e b,a - Here, R b,a is a phase, but in general, R b ' a is a unitary matrix. 
We should mention that in general R b,a is not the inverse of Rc ,b . Their product 
involves the twists of particles. 



As we have seen before anyons can be fused or splitted, therefore braidings should 
be compatible with them. For example, given two anyons c, d, we may first split d 
to a, b, then braid c with a followed by braid c with b, or we may braid c and d first, 
then split d into a, b. These two processes are physically equivalent, therefore their 
resulting matrices should be the same. Applying the two operators on the fusion 
tree basis e%f, we have an identity in pictures: 



b 




a 



b 



c 



a 



b 



c 




d 
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Fig. 4. (color online) The hexagon relation for "i?-moves" and "f-moves" . 

The same identity can be also obtained as a composition of F-moves and braid- 
ings as shown in Fig. |U It follows the composition of the 6 matrices, hence the 
name hexagon, should be the same as the identity. The resulting equations are 
called hexagons. There is another family of hexagons obtained by replacing all right- 
handed braids with left-handed ones. In general, these two families of hexagons are 
independent of each other. Similar to the pentagons, a mathematical theorem says 
that the hexagons imply all other consistency equations for braidings. 

Written explicitly for Fibonacci anyons the hexagon equation reads 



R T C ' T (Fl 



Rl, 



(Frr b K' b (Frr a 



b 



(2-5) 



where again the indices o, b, c label the internal edges of the fusion trees as shown 
in Fig. HJ Inserting the .F-matrix (|2-4[) and realizing that braiding a particle around 



the trivial one is trivial: Rr' 1 = Rt' t = 1 the hexagon equation becomes 



RYRV^>- xt2 



RYR T r T ^- l/2 

-(^)V 1 



i#V-i + ^ 2 (i - 
(i - RV)v~ V2 R T,T 



R T r 



7T V~ 3/2 



which has the solution 



3 +47ri/5 



Rl 



-3ni/5 



(2-6) 



(2-7) 



The combined operation of a basis transformation F before applying the ii-matrix 
is often denoted by the braid-matrix B 



B = F? TT R TT F t 



arr 
c 



(2-t 
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Using a basis {|a&c)} for the labelings adjacent to the two anyons to be braided the 
basis before the basis transformation is 

{|lrl), \tt1), \1tt), \t1t), \ttt)} (2-9) 

and after the basis change to a basis {|afec)} using an F matrix the basis is 

{(111), \tt1),\1tt), \t1t), \ttt)}. (2-10) 

as illustrated here. 



T T F 



\b lb 



a b c 



a c 



a c 



a b 



In this representation the F-matrix is given by 

/ 1 



F 



\ 



1 



1 



\ 



-V2 _ M -1 



^ ) 



and the i?-matrix is 

R = diag(e 4 ^ /5 , e~ 3ni/5 , e~ 3wi/5 , e^.e" 3 "/ 5 ) 
We finally obtain for the braid matrix 

/ e 47ri/5 



B = FRF' 



-37ri/5 



-3-7ri/5 



(2-11) 



(2-12) 



\ 



^-l e - 47 ri/5 _^-l/2 e -27ri/5 



V 



_^-l/2 e -27ri/5 



-if 



J 

(2-13) 

With the explicit matrix representations of the basis transformation F and the 
braid matrix B, we are now fully equipped to derive matrix representations of Hamil- 
tonians describing interactions between Fibonacci anyons. 



§3. Hamiltonians 



Considering a set of Fibonacci anyons we will now address how to model inter- 
actions between these anyons. Without any interactions, the collective state of the 
anyons will simply be described by the large degenerate manifold described by the 
Hilbert space introduced above. If, however, the anyons interact, this degeneracy 
will be split and a non-trivial collective ground state is formed. In this section we 
will first motivate a particular type of interaction, which generalizes the well-known 
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Heisenberg exchange interaction to anyonic degrees of freedom, and then explic- 
itly derive various Hamiltonians of interacting Fibonacci anyons that correspond to 
well-known models of SU(2) spin-l/2's. 

Two SU(2) spin-l/2's can be combined to form a total spin singlet or a total 
spin triplet 1, which in analogy to the anyonic fusion rules we might write as 

1/2 1/2 = 00 1. 

If the two spins are far apart and interact only weakly, these two states are degen- 
erate. However, if we bring the two spins close together a strong exchange will be 
mediated by a virtual tunneling process and the degeneracy between the two total 
spin states will be lifted. This physics is captured by the Heisenberg Hamiltonian 
which for SU(2) spins is given by 

"^Heileliberg = ' &3 = ~^ ~ ^ ~ ^f) = "2 I X/ ^ij ~ 2 j ' ^ 

<*?'} (ij) ) 

where S% and Sj are SU(2) spin 1/2's, T^- = Si + Sj is the total spin formed by the 
two spins Si and Sj, a (uniform) coupling constant is denoted as J, and the sum 
runs over all pairs of spins i, j (or might be restricted to nearest neighbors on a given 
lattice). Of course, the Heisenberg Hamiltonian is just a sum of projectors IJ^ onto 
the pairwise spin singlet state as can be easily seen by rewriting the spin exchange 
Si ■ Sj in terms of the total spin Tij. Antiferromagnetic coupling (J > 0) favors an 
overall singlet state {Tf- = 0), while a ferromagnet coupling (J < 0) favors the triplet 
state (ig = 2). 

In analogy, we can consider two Fibonacci anyons. If the two anyons are far 
apart and weakly or non-interacting, then the two states that can be formed by 
fusing the two anyons will be degenerate. If, however, the anyons interact more 
strongly, then it is natural to assume that the two fusion outcomes will no longer 
be degenerate and one of them is energetically favored. We can thus generalize the 
Heisenberg Hamiltonian to anyonic degrees of freedom by expressing it as a sum of 
projectors onto a given fusion outcome 

-T/Fibonacci j \ rjl /o n\ 

rl Heisenberg — J 2^t V V J ' Z J 

where U\j is a projector onto the trivial channel. 
3.1. The golden chain 

We will now explicitly derive the matrix representations for simple models of 
interacting Fibonacci anyons. In the simplest model we consider a chain of Fibonacci 
anyons with nearest neighbor Heisenberg interactions as shown in Fig. [6j 
This Hamiltonian favors neighboring anyons to fuse into the trivial (1) channel by 
assigning an energy — J to that fusion outcome. To derive the matrix representation 
in the fusion tree basis we first need to perform a basis change using the F-matrix 
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SU(2) spins 



Weak interaction 



\ 1 /2 1/2 = 1 / 
1/2 1/2 = / 

Strong interaction 

1/2 1/2 = 1 



A oc J 



1/2 1/2 = 



Fibonacci anyons 



Weak interaction 



T 



T T = T 



T <g> T = 1 



T 



Strong interaction 
r r = r 



T T 



A oc J 
= 1 



Fig. 5. (color online) The generalized Heisenberg model: While for two weakly interacting SU(2) 
spin-l/2's (left panel) the total singlet and triplet states are degenerate, strong interactions will 
lift this degeneracy. The Heisenberg Hamiltonian explicitly opens a gap of the order A oc J be- 
tween the two states. Applying a similar idea to the effect of interactions between two Fibonacci 
anyons (right panel), a generalized Heisenberg Hamiltonian will energetically distinguish the two 
fusion outcomes. 



T T T T T T 



abed 



Fig. 6. (color online) The "golden chain" of pairwise interacting Fibonacci anyons. The interactions 
are indicated by ellipses around the anyons. Our choice of fusion path basis is indicated in the 
right panel. 



T T 



a 



c 




T T 

a 



and then the Hamiltonian is just FII F with the projector 77 given by 77 = 
diag(l,0, 0, 1,0) in the basis f|2- 1Q|1 . Written explicitly this becomes 



-J{FD\{FDl 



(3-3) 
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The matrix representation in the basis (I2-9P then reads 



/ 1 



\ 







- J 



-3/2 
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3.2. Three-anyon fusion 

For the second model we include longer range interactions, preferring now to 
fuse three adjacent anyons into the trivial particle. For this we have to perform two 
basis changes to obtain the total charge of three anyons as shown here: 



T T T 




T T T 

b4= 




abed 



a 



d 



a d 



The basis states of three anyons are given by the labelings of the four edges 
\abcd) between and adjacent to the three anyons: 



{|1tt1), |1t1t), \1ttt), |rlrl), \ttt1), \tt1t), tItt), \tttt)} . 
The first matrix is 



/ 1 



Fi 



1 



-1/2 



-1/2 



<P 



\ 



-1/2 



r 



-1/2 
^ ) 



changing to the basis \abcd) 

{jlrrl}, |lllr), \1ttt), \t1t1), \ttt1), \tt1t), tItt), \tttt)} 
and then a second basis change with 

/ 1 



\ 



1 



1 



1 



<p~W 
1 
tp- 1 ' 2 ) 



to a basis abedj 



(3-5) 



(3-6) 



(3-7) 



(3-8) 



{ | 1t1 1) , |11tt), \1ttt), |t1t1), \ttt1), \t1tt), ttIt), \ tttt)}. (3'9) 
Combined with the projector P3 = diag(l, 0, 0, 0, 0, 1, 0, 0) the Hamiltonian matrix 
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then becomes 



/ 1 



-J3F1F2P3F2F1 — — J3 



V 





-2 


¥?" 


-2 


-<^/ 2 


¥>~ 


-2 




-2 




-ip- 


-5/2 




-5/2 





/ 

(3-10) 



3.3. Next-nearest neighbor interactions 

Another possibility for longer-range interactions is to fuse two particles at larger 
distance. Here we consider the case of next-nearest neighbor interactions. To deter- 
mine the fusion result of these two anyons we first need to bring them to adjacent 
positions by braiding one of the particles with its neighbor as shown here: 





r r r 



In contrast to Abelian particles, such as ordinary spins, for non-Abelian anyons 
it matters whether we braid the two anyons clockwise or counter-clockwise, since 
the braid matrix (for braiding clockwise) is different from its inverse (braiding anti- 
clockwise). Indicating the interaction of two anyons by loops around the two, the 
two ways of braiding correspond to the anyons fusing above or below the one between 
then as shown here: 




Using the same basis as above we obtain for the Hamiltonian for upper and lower 



interactions 

#2,abovc = BH X B~ X = H = FPF 



( 



J 2 ^ 3/2 



<P 



-1/2 

„-2 



<P 



r 

1/2 



-r 



1/2 
2 



r 

1/2 



_ (/9 -l/2 r -4 



Lp~^r 2 



(3-11 



<P 



-l r 2 



— 1 —2 

tp L r 



-3/2 
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and 

#2,below = B~ 1 HiB = #2,above ) (3T2) 

where B is the braid matrix acting on the first two anyone (first three labels) and 
r = exp(i7r/5). For the symmetrized model containing both terms we have 

(3-13) 
\ 

2^-2 _p-5/2 

299-2 99- 7 /2 
-if- 1 2if~ 2 if" 7 / 2 
^-7/2 2 <p- 3 J 

As for SU(2) spin chains, the three-particle interaction can be written as a sum 
of nearest and next-nearest neighbor terms: For SU(2) spin-S this is 

#3 = (Si + S 2 + S 3 ) 2 = 3S(S + 1)1 + 2S 1 S 2 + 25i5 3 + 2S 2 S 3 

= 3S(S + 1)1 + 2H\ 2 + 2Hf + 2H 2 13 (3-14) 

where the H± indicates a nearest neighbor interaction between sites i and j and H l 2 ° 
a next-nearest neighbor one. Mapping an ^3 chain to a chain containing both Hi 
and H 2 terms we find J\ = 4J3 and J 2 = 2J3. The H 3 chain thus corresponds to the 
Majumdar-Ghosh chain 9 ) with an exact ground state of singlet dimers. 
Similarly we find for Fibonacci anyons 

H 3 = -1 - ip- 2 H{ 2 - <p- 2 Hf 3 - (/? _1 #2 13 (3-15) 

and thus J\ = — 2ip~ 2 J 3 and J 2 = — ip^ 1 J3. Again the pure H 3 chain corresponds 
to the Majumdar-Ghosh chain with an exact ground state of anyonic dimers fusing 
into the trivial particle. 10 ) 

Generalizations of the next nearest neighbor interaction to longer distances is 
straightforward, the only complication arises from the various possible orientations 
of the braids, giving 2 r_1 different terms at distance r. 

3.4. The two-leg ladder 

The final model we will present is a two-leg ladder consisting of two coupled 
chains. Unlike the case of the chain, where it was natural to just use the standard 
fusion tree as basis, there are several natural choices here. Choosing the zig-zag 
fusion path indicated by a line in Fig. [7] minimizes the interaction range on the 
fusion tree. 

Distinguishing rung interactions with coupling constant Jj_ and chain interac- 
tions with coupling J, we see that the rung terms are just nearest neighbor interac- 
tions and the chain terms are next nearest neighbor interactions above or below the 



H 2 — i?2,below + -f^2,above 

/ 

2<p~ 2 -^l 2 

-ip- 5 / 2 2<p~ l 

-J2 
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ft 





Fig. 7. (color online) A two-leg ladder model, consisting of two chains of Fibonacci anyons. Pairwise 
interactions are indicated by ellipsis around the anyons. Our choice of fusion path is indicated 
by a line in right panel. 



other anyon for the upper and lower chain respectively: 

iw = j±E ^ + J E + J E ^Ste • (3-i6) 

i i i 

Having established the explicit matrix representations of these Hamiltonians we 
can now analyze their ground states thereby addressing the original question of what 
kind of collective ground states a set of Fibonacci anyons will form in the presence 
of these interactions. 



§4. Alternative SU(2)fc formulation 



Before turning to the collective ground states of the Hamiltonians introduced 
above, we will describe an alternative reformulation of the "golden chain", Eq. (j3- 2D . 
which allows us to generalize the latter and, in particular, allows us to view the 
golden chain and its variants as certain deformations of ordinary SU(2) spin chains. 
Technically, this generalization is based on connections with the famous work by V. 
F. R. Jones 13 ) on representations of the Temperley-Lieb algebra. 14 ) 

Specifically, this reformulation is based on the fusion rules of the anyon theory 
known by the name 'SU(2) at level k\ which in symbols we denote as SU(2)fc, for 
k = 3. For arbitrary integer values of the parameter k there exist particles labeled 
by 'angular momenta' j = 0, 1/2, 1, 3/2, k/2. The result of fusing two particles 
with 'angular momenta' j\ and ji yields particles with 'angular momenta' j where 
j = \ji~h\, \ji-h\ + l, ■■■,mm{j 1 +j 2 ,k-(j 1 +j 2 )}, eachoccuring with multiplicity 
^jij2 = ^' When k — > oo this represents ordinary SU(2) spins, for which all possible 
values of (ordinary) SU(2) angular momenta j = 0,1/2,1,3/2,.... appear, and the 
fusion rule turns into the ordinary angular momentum coupling rule. Finite values 
of k represent a 'quantization' of SU(2), amounting to the indicated truncation of 
the range of 'angular momenta'. Thus, our reformulation will allow us to view the 
"golden chain" as a deformation of the ordinary SU(2) Heisenberg spin chain by a 
parameter 1/k. For the special case where 1/k = 1/3 we obtain, as we will now 
briefly review, the chain of Fibonacci anyons introduced in Eq. f|3-2|) above. 

Our reformulation of the "golden chain" is based on two simple observations: 
(i) first, one recalls that at k = 3 (where the four particles j = 0, 1/2, 1,3/2 exist) 
the fusion rule of the particle with j = lisl<8>l = 0©l (all entries label 'angular 
momenta'), which is identical to that of the Fibonacci anyon r, e.g. r <g> r = 1 © r. 
Note that the trivial particle, previously denoted by 1, is now denoted by 'angular 
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momentum' j = 0. (ii) Secondly, we recall that the particle j = 3/2 can be used 
to map the particles one-to-one into each other (it represents what is known as an 
automorphism of the fusion algebra), i.e. 

3/200 = 3/2, 3/2 1/2 = 1, 3/20 1 = 1/2, 3/2 3/2 = 0. (4-1) 

We denote this automorphism by 

3/20j=j: j^j. (4-2) 

Note that when j is an integer, then j is a half-integer, and vice-versa. 

Using the observation (i) above we may first write an arbitrary state vector in 
the Hilbert space which can be represented pictorially as 

T T T T T T 



T 1, r 1, r 1, t 



in the new notation using the 'angular momenta' as follows 

111111 



1 0,1 0,1 0,1 .. . 

Let us now consider the matrix element of the projector [discussed in Eq.s (|3-2|) 
-(USD above] 

(...a',b',c...\II 1 \...a,b,c...) (4-3) 

where ...a' ,b' , c' , ...a,b,c, ... G {j = 0,j = 1}. A moment's reflection shows that in 
fact a' = a and c' = c. This matrix element is graphically depicted in the leftmost 



b' 



1 




i 




n 1 




1 




i 



a b c 
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1 

1/2 



1 



_y_ 

3/2 



"3/2" 



1 

1/2 
1 



n 1 



i 

1/2 ' 
1 



3/2 



"3/2" 



1 

1/2 
1 



a V 
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n 1 




1 
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2 
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(I 




b 




C 



a 



Fig. 8. (color online) Transformation of the projector II 1 . 
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picture in Fig. [8j where the direction from the right (ket) to the left (bra) in (14-3|) is 
now drawn vertically upwards from bottom to top in the Figure. 

We may now perform a number of elementary steps using a 1 x 1-dimensional F- 
matrix (which is just a number) involving the Abelian particle j = 3/2. Specifically, 
we perform the steps depicted as 



11 11 




a b c a b c a b c 



on both, the upper and the lower legs of the projector in Fig. [8l where in the last 
step we use 

3/2 3/2 3/2 3/2 

3/2 ,■ 3/2 3/2 ,. 3/2 Y' ^ 

3/2 3/2 3/2 3/2 | / \ 

3/2 3/2 3/2 3/2 



Subsequently we use, in the diagram in the middle of the projector equation Fig. 
the transformation 

1/2 1/2 1/2 1/2 1/2 1/2 




where in the first step we have applied the relation 



1 3 

2 2 1/2 3/2 




a b abb 



The projector 77 1 with four j = 1/2 particle legs (depicted in the rightmost diagram 
of Fig. [8]), which we obtain at the end of this sequence of elementary steps is the 
central object of interest to us. 
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As it is evident from the rightmost picture in Fig. [8l the projector 77 1 acts on 
basis states of the form 

...a,b,c,d,e,...^ , (4-4) 

that is on states in which integer and half-integer labels are strictly alternating. Such 
basis states can then be depicted as 

111111 

2 2 2 2 2 2 



1 I ? o,l I 3 0,1 .. . 

2' 2 2' 2 

Fig. 9. 'Angular momentum' description of the wavefunction for the Fibonacci theory SU(2)3. 

The projector 77 1 that we have thereby obtained is a central and well known 
object of study in Mathematics: Specifically, let us form the operator 

ei := d n\ (4-5) 

where d = <p = (y/b + l)/2 is the golden ratio. The index i indicates the label along 
the fusion chain basis of 'angular momenta' in (|4-4p on which it acts. For example, 
the operator 77 1 in (|4-3p as drawn in Fig. [8] acts on the label b which is, say, in the i th 
position in the chain of symbols characterizing the state. With these notations, the 
operators e% are known to provide a representation of the Temperley-Lieb algebra 14 ) 

6j — ds-i , 

[ei,ej] =0, \i — j\ > 2. 

with d-isotopy parameter d = ip = (v5 + l)/2. The so-obtained representation of the 
Temperley-Lieb algebra is only one in an infinite sequence with different values of the 
d-isotopy parameter, discovered by V. F. R. Jones. 13 ) Specifically, in our notations, 
this infinite sequence is labeled by the integer k, specifying the (i-isotopy parameter 
d = 2cos[-7r/(fe + 2)]. For any integer value of k, the Jones representation is given 
by the matrix elements of the operator in the basis of type (|4-4|) . Recall from the 
discussion at the beginning of this paragraph that for general values of the integer 
k the labels ...a, b, c... are the 'angular momenta' in the set j = 0, 1/2, 1, 3/2, k/2. 
Thus it is in general more convenient to write these basis states as 

\- ■ -Ji-iJiJi+i ■ ■ ■) (4-6) 

with ji-i, ji, ji+i G {0, 1/2, 1, 3/2, k/2}; the sequence must satisfy the condition 
that ji is contained in the fusion product of with an 'angular momentum' 1/2, 
ji+i is contained in the fusion product of ji with 'angular momentum' 1/2, etc. - 
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this generalizes the states depicted in Fig. [9] (which is written for k = 3). Within 
this notation the matrix elements of the operator are (compare I4-3P 



(■••ii-i)ii)ii+i---l e il---ii-i)iijii+i---) — ^-1,^+1 



5050 

\ qO qO ' 
(4-7) 



where 



2 . (2j + l)(2/ + l) 1 

5 ; = V(^) sm[7r fcT2 ] " (4 ' 8) 



This, then, defines a generalization of the original "golden chain", for which 
k = 3, to arbitrary values of k. It is interesting to note that in the limit k — > oo, the 
so-defined generalized anyonic chain turns precisely into the ordinary SU(2) spin- 
1/2 Heisenberg chain. Therefore, these generalized "golden chains" (for arbirary 
integer values of k) provide, as alluded to in the introduction of this section, a 
certain generalization of the ordinary SU(2) spin-1/2 Heisenberg chain. Therefore, 
it is natural to ask questions about the behavior of ordinary spin-1/2 chains, in 
the context of these generalized golden chains. Indeed, as we have described in 
various publications 10 ^ 12 -* and as reviewed in this note, many of the physics questions 
familiar from the ordinary Heisenberg chain find their mirror-image in the behavior 
of our generalized golden chains. The following section is intended to give a brief 
flavor of these results and parallels. 

Before we proceed, let us briefly mention that for 'antiferromagnetic' coupling 
of the generalized spins of the golden chain (for general k), meaning that projection 
onto the singlet state (trivial anyon fusion channel) between two spins is energetically 
favored, one obtains a gapless theory which turns out to be precisely the (k — l)-th 
minimal model of conformal field theory 15 ) of central charge c = 1 — 6/[(k+l)(k + 2)]. 
(For k = 2 this is the Ising model, for k = 3 the tricritical Ising model, and so on.) 
In the opposite, 'ferromagnetic' case, meaning the case where the projection onto 
the 'generalized triplet' state of two neighboring generalized spins is energetically 
favored, one obtains another well known sequence of gapless models: these are the 
Zfr parafermionic conformal field theories, 16 ) of central charge c = 2(k — l)/(k + 2), 
which have more recently attracted attention in quantum Hall physics as potential 
candidates for certain non-Abelian quantum Hall states, known as the Read-Rezayi 
states. 5 -* For a summary see Table I. 



§5. Collective ground states 



In this final section we will round off the manuscript by shortly reviewing some 
recent numerical and analytical work analyzing the collective ground states formed 
by a set of Fibonacci anyons in the presence of the generalized Heisenberg interactions 
introduced in the previous section. 

For the "golden chain" investigated in Ref. 11) , a one-dimensional arrangement 
of Fibonacci anyons with nearest neighbor interaction terms, it has been shown (as 
already mentioned above) that the system is gapless - independent of which fusion 
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level k 


AFM 


FM 


2 


c= 1/2 


c= 1/2 




Ising 


Ising 


3 


c = 7/10 


c = 4/5 




tricritical Ising 


3-state Potts 


4 


c = 4/5 


c = 1 




tetracritical Ising 


Zzi-parafermions 


5 


c = 6/7 


c = 8/7 




pentacritical Ising 


Zs-parafermions 


k 


c=l-6/[(fc + l)(fc + 2)] 


c = 2(fc-l)/(fc + 2) 




fc t,l -multicritical Ising 


Z fc -parafermions 


oo 


c = 1 


c= 2 


Heisenberg AFM 


Heisenberg FM 



Table I. The gapless (conformal field) theories describing the generalized spin-1/2 chain of anyons 
with SU(2) fe non-Abelian statistics. 

channel, the trivial channel (1) for antiferromagnetic coupling, or the r-channel for 
ferromagnetic coupling, is energetically favored by the pairwise fusion. The finite- 
size gap A(L) for a system with L Fibonacci anyons vanishes as A(L) oc (1/L) Z=1 
with dynamical critical exponent z = 1, indicative of a conformally invariant energy 
spectrum. The two-dimensional conformal field theories describing the system have 
central charge c = 7/10 for antiferromagnetic interactions and c = 4/5 for ferro- 
magnetic interactions, respectively, corresponding to the entry for k = 3 in Table 
I. In fact, a direct connection to the corresponding two-dimensional classical mod- 
els, the tricritical Ising model for c = 7/10 and the three-state Potts model for 
c = 4/5, has been made: Realizing that (as reviewed briefly in the previous section) 
the non-commuting local operators of the "golden chain" Hamiltonian form a well 
known representation 13 ) of the Temperley-Lieb algebra 14 ) (with <i-isotopy parameter 
d = <p), it has been shown 11 ) that the Hamiltonian of this quantum chain corresponds 
precisely to (a strongly anisotropic version of) the transfer matrix of the integrable 
restricted-solid-on-solid (RSOS) lattice model, 17 ) thereby mapping the anyonic quan- 
tum chain exactly onto the tricritical Ising and three-state Potts critical points of 
the generalized hard hexagon model. 11 )' 18 ) A corresponding exact relationship holds 
in fact true for the chains at any value of the integer k (Table I) . 

While this correspondence of the "golden chain" with the special critical points 
of the classical models might initially seem accidental, it turns out that the quantum 
system exhibits in fact an additional topological symmetry that actually stabilizes the 
gaplessness of the quantum system, and protects it from local perturbations which 
would generate a gap. In particular, it was shown that all translational-invariant 
relevant operators that appear in the quantum system, e.g. the two thermal operators 
of the tricritical Ising model (antiferromagnetic case) e, e' with scaling dimensions 
1/5 and 6/5 respectively, are forbidden by this topological symmetry. 11 ) 

A more detailed connection to the underlying two-dimensional classical models 
has been made by considering the effect of a competing next-nearest neighbor in- 
teraction or, equivalently, a three-anyon fusion term in the anyonic analog of the 
Majumdar-Ghosh chain, 10 ) derived explicitly above. The rich phase diagram of this 
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exact 
solution 




tan = <j>/2 

critical 3 -state Potts 
\c=4/5 with S 3 -symmetry 



te 



exact 
solution 



0=3ji/2 tetracritical Ising 
Majumdar-Ghosh c=4/5 



Fig. 10. (color online) The phase diagram of the anyonic Majumdar-Ghosh chain, for Fibonacci 
anyons, as presented in Ref. 10). The exchange couplings are parametrized on the circle by 
an angle 6, with a pairwise fusion term J2 = cos# and a three-particle fusion term J3 = sm9. 
Besides extended critical phases around the exactly solvable points (9 = 0, 71") that can be 
mapped to the tricritical Ising model and the 3-state Potts model, there are two gapped phases 
(grey filled). The phase transitions (red circles) out of the tricritical Ising phase exhibit higher 
symmetries and are both described by CFTs with central charge 4/5. In the gapped phases 
exact ground states are known at the positions marked by the stars. In the lower left quadrant 
a small sliver of an incommensurate phase occurs and a phase which has ^-symmetry. These 
latter two phases also appear to be critical. 



model is reproduced in Fig. [TOj Besides extended critical phases around the "golden 
chain" limits (9 = 0, ir) for which an exact solution is known, there are two gapped 
phases (grey filled) with distinct ground-state degeneracies and well-defined quasipar- 
ticle excitations in the spectrum. The tricritical Ising phase ends in higher-symmetry 
critical endpoints, with an S^-symmetry at 9 = 0.1767T and an Ising tetracritical point 
at 9 = — 0.472-7T. The transitions at these points spontaneously break the topological 
symmetry and in the case of the tetracritical point also the translational symmetry 
thereby giving rise to two- fold and four- fold degenerate ground states in the adjacent 
gapped phases, respectively. For an in-depth discussion of this phase diagram we 
refer to Ref. 10). 

Finally, the effect of random interactions on chains of (Fibonacci) anyons has 
been studied in Refs. 19), 20). For random, 'antiferromagnetic' interactions the 
random system is found to flow to strong disorder and the infinite randomness fixed 
point is described by a generalized random singlet phase. 19 ) For a finite density of 
'ferromagnetic' interactions an additional 'mixed phase' infinite randomness fixed 
point is found. 20 ) 
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